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Definition of the phase of oscillations is straightforward for deterministic periodic processes but 
nontrivial for stochastic ones. Recently, Thomas and Lindner in [Pliys. Rev. Lett., v. 113, 
254101 (2014)] suggested to use the argument of the complex eigenfunction of the backward density 
evolution operator with the smallest real part of the eigenvalue, as an asymptotic phase of stochastic 
oscillations. Here I show that this definition does not generally provide a correct asymptotic phase. 

PACS numbers: 05.40.-a 


Notion of the phase of periodic oscillators lies at the 
heart of characterization of oscillatory processes. In de¬ 
terministic systems, the phase on the limit cycle is de¬ 
fined straightforwardly as a 27r-periodic variable on the 
cycle which grows in time uniformly. For stable os¬ 
cillations, the definition can be extended also to the 
basin of attraction of periodic regime, via construction 
of isochrons - the surfaces of Poincare sections of the 
flow with the return time being exactly the period of 
oscillations. The isochrons are surfaces of the constant 
phase (sometimes called also the asymptotic phase). Re¬ 
cently, in [lj the notion of isochrons has been extended to 
stochastic systems, where they have been defined as the 
surfaces of constant mean first return time. In paper [2j, 
a different definition of the asymptotic phase of stochastic 
oscillations has been suggested, based on the properties 
of the operator describing the evolution of the probabil¬ 
ity density (we assume for simplicity of presentation that 
it is the Fokker-Planck operator). If the nontrivial eigen¬ 
value of the backward Fokker-Planck operator with least 
negative real part is complex, the systems has been called 
robustly oscillatory. With some further minor technical 
conditions, for such systems the asymptotic phase was de¬ 
fined in Q as the argument of the complex eigenfunction 
corresponding to the first nontrivial complex eigenvalue. 

In this comment I present a simple analytically solv¬ 
able example, where the definition of [2] does not pro¬ 
vide a proper phase. I consider a limit cycle in a three- 
dimensional phase space, and use poloidal coordinates 
(9,cf>,r) around it. The limit cycle corresponds to r = 0 
and is parametrized by its phase 9 , which is governed 
by 9 = f2, where f2 is the frequency of the oscillations. 
Transversal to the limit cycle direction is spanned by vari¬ 
ables (r, (f>). The limit cycle is stable, and the transversal 
perturbations decay according to f = —Sr, <f> = ui, where 
—S is the rate of decay of perturbations. In the whole 
vicinity of the limit cycle, the variable 9 is the proper 
phase. 

We assume now that these oscillations are subject to 


independent Gaussian white noise terms in all variables: 

9 = fl + crg^g(t) , 

<j> = uj + cr^it) , 
f = -Sr + oy£ r (f) ) 


where (^e(t)^e(t')) = 2 S(t — t') and similarly for other 
terms. If the noise terms are sufficiently small, the oscil¬ 
lations are slightly perturbed deterministic ones, and are 
well defined for all relations between the noise intensi¬ 
ties; the approach of Ref. Q yields here the proper phase 
9 = const. 

Let us apply approach of Ref. Q. An eigenvalue prob¬ 
lem for the backward Fokker-Planck equation 
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can be solved by separation of variables. The eigenfunc¬ 
tions 


Pnmi = e ine Ri(r) 

are parametrized by three integers ( n,m,l ), and the 
eigenvalues are 

A nmi = inCl+imw—nog—mOf—pi, n,m,l = 0,1,2,... , 
where /3; are the eigenvalues of the radial equation 

<J r~^ L ~ r ^ Rl + & lRl = 0 ' 

The latter is nothing else but the equation for the 
Ornstein-Uhlenbeck process, with the eigenfunctions be¬ 
ing Hermite polynomials and the eigenvalues fii = IS (see 
0, section 5.2.6 for details). 

According to [§], we focus on the nontrivial eigenvalue 
with least negative real part. First, we observe that for <5 
small compared to , cr 2 , this eigenvalue is A 001 and real, 
so that the process is not robustly oscillatory according 
to classification of Ref. Q. If parameter <5 is large, then 
this eigenvalue is complex, and one has two possibilities: 
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1. If <T 0 < ere, then X min = A 0 io = ico - a^. In this 
case the eigenfunction is poio = Ro(r)e 1 ^. 

2. If <j 0 > ag, then A m i„ = Aioo = ifl — <jg. In this 
case the eigenfunction is p W o = Ro{r)e ie . 

One can see that only in case (2) the proper phase 9 
is recovered, while in case ( 1 ) the variable </>, which is 
transversal to the correct phase 9 1 is delivered as an 
asymptotic phase variable according to the method of 
Ref. 0. 

In conclusion, while the phase is a rotating variable, in 
the phase space of a noisy dynamical system there can be 
many rotations corresponding to complex eigenvalues of 


the density evolution operator. Therefore one generally 
cannot identify the phase using the least stable complex 
eigenfunction. 
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